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ABSTRACT 

Dye densities of multi-layered films are determined by applying a 
regression analysis to the spectral response of the composite transparency 
The amount of dye in each layer is determined by fitting the sum of the 
individual dye layer densities to the measured dye densities. From this, 
dye content constants are calculated. 

Methods of calculating "equivalent exposures" are discussed. 
"Equivalent exposures" are a constant amount of energy over a limited 
band-width that will give the same dye content constants as the real 
incident energy. Methods of using these equivalent exposures for analysis 
of photographic data are presented. 


INTRODUCTION 


All the information recorded at each point on multi-layered films 
can be characterized by one constant for each layer of the film. The 
constant for each layer is the number a normalized dye density must be 
multiplied by in order to equal the dye in that layer for a particular 
scene response. 

All multi-layered films have characteristic colors associated with 
each layer of the film. The spectral properties of the dyes that make up 


these colors are well known^, though not necessarily stable with time. 
An example of the spectral characteristics for Kodak Ektachrome MS 
Aerographic film type 2448 is given in Figure 1. In this case the film 
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has three layers as do most multi-layered films* The colors formed in 
the three layers are yellow, magenta and cyan. The properties of this 
film are well known. 2 The method of analysis is developed using type 
2448 as an example though the methods can be easily generalized to any .. 
multi-layered film. 

The spectral response or color at each point on a transparency is 
due to the dyes in the three layers at that point. The spectral properties 

t 

of the dyes in each of the layers at a point can be represented by a 
constant (one for each dye) times a normalized dye density at each wave- 
length. Figure 1 shows the three dye density curves for film type 2448. 

Thus the spectral properties of the transmitted light through the developed 

film transparency are determined by a linear superposition of constants 

times these dye density curves. Mathematically this can he written. 

D(x,y,X) = Ki(x,y)Y(A)+ K 2 (x,y)M(A) + K 3 (x,y)C(A) + B(A) (1) 

and I (x , y , A ) = I q (x,y ,X) 10 ” :D Cx»y,X) (2) 

where D(x,y,X) is the density at any point as a function of x, y and 
of wavelength (X) 

Y(A) is the normalized dye density for the yellow layer as a 
function of wavelength 

M(A) is the normalized dye density for the magenta layer as a 
function of wavelength 

C(A) is the normalized dye density for the cyan layer as a 
function of wavelength 

B(A) is the density of the film base as a function of wavelength 

Ki(x,y), K 2 (x,y) » K 3 (x,y) are constants related to the amount of 
dye on the three respective film layers and only a function of 
position on the film format (dye content constants) 
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I(x,y,X) is the transmitted light through the transparency at 

each point on the film and is a function of wavelength 

l 0 (x,y,*) is the incident light on the transparency 

At any given point, the quantities experimentally determined are ^ .. 

usually the total transmitted light and the incident light. From this 

the density can be calculated for that point. Because these quantities 

are a function of the wavelength of the light (color) , they must be 

measured at a single wavelength. Actually, for this treatment, measurements 

must be acquired for at least three different wavelengths at the point of 

interest on the transparency. Three equations can be written and solved 

for the three dye content constants. It is then possible to determine an 
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equivalent integral spectral density or analytical spectral density at 
any wavelength. 

The transmittance of a typical film base is represented in Figure 2. 
Predicted and observed values of transmittance for colored objects are 
shown in Figures 3, 4 and 5. 

Equivalent 'Energies 

An equivalent energy is defined as a constant amount of energy through 

out the spectrum that will give the same dye constant in a film layer as the 

real incoming radiation. There will be a separate equivalent energy for 

each layer of the film. These equivalent energies can be calculated by 

using the known sensitivities of the film layers. 

t 

Let: S i (X) be the sensitivity of the i layer of the film as 
a function of wavelength 

2 

e (A) be the radiation falling at a point on the film (ergs/cm / ) 

th 

E. be the equivalent energy for the i film layer 
i = 1 for the yellow layer 
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Figure 5. Transmittance vs. wavelength of a red object. 
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i = 2 for the magenta layer 


i = 3 for the cyan layer 
Mathematically we can write: 

/E.S.(X)dX - /e(X)S, (X)dX - (3) 

IX X 

where the integral is over all wavelengths. 

The units of S_, (X) are in l/(ergs/cm ) to give a density of 1 with 
i 

standard processing. Thus the integral of radiation times the sensitivity 
will give the fraction of exposure required to give a final density of 1 


with standard processing for that radiation. 

Since E on the left of (3) is a constant we can write (3) as: 
i 

/e (A)S-j (A)dA 
i / S.(A)dA 


(4) 


In essence, this E^ is a weighted average of incoming radiation over the 
sensitivity of the layer. 

In themselves these E's are not very useful. They represent three 
different constant energy spectrums. More useful would be one equivalent 
energy spectrum that can replace the incoming radiation spectrum and give 
the same dye densities. Using the same type of treatment, this equivalent 
energy spectrum can be found in terms of the E’s and thus in terms of 
incident energy spectrum s(A) and the sensitivity of the film layer, S^(A) . 

For color and color IR films these three equivalent energies are 


designated as EB, EG and ER and EG, ER and El respectively. An example 
of these equivalent energies is graphically represented in Figure 3. EB* 
is a constant amount of energy between .4 and .5 microns, EG is a constant 
amount between .5 and .6 microns, ER is a constant amount between .6 and 
.7 microns, and El is a constant amount between .7 and .9 microns. 

The actual values of the three equivalent energies for each of the 
film types can be easily found by substituting the equivalent energy spectrum 

*Kodak filter 2B must be used for it if EB is to correspond to the range 
.4 to ,5p. 



into equation (4) in the place of e(A), being careful to integrate over 


the proper limits. In the case of color film this reduces to: 

/; 5 EB S.(X)dA + /; 6 EG S . (A)dA + /; 7 ER S. (A)dX 
. 4 1 *5 i * p 1 

;; 7 s, (x)dx 

.4 i 


E. » 

l 


( 5 ) 


Equation (5) is really 3 equations that can be v/ritten explicitly as: 

( 6 ) 

(7) 


E 1 A il EB + A 12 EG + A 13 ER 


E 2 = A^EB + A^EG + A^ER 
E 3 - A^EB + A' 2 EG + A« 3 ER 


( 8 ) 


where 


A , . /; 5 s 1 (x)dx 

A 11 A— 

/^ / S 1 (X)dA 


. 6 , 


., /: w s 1 (x)dx 

A ' = . 5 1 


A 1 

21 


A’ ■ 
31 


/^ 5 S 2 (X)dX 

/^ 7 S 2 (X)dX 

/• 5 S 3 (X)dX 

/• 7 S 3 (X)dX 


12 


A’ 

A 22 


A' = 
A 32 


. {e 7 s 1 (x)dx 


;^ 7 s 1 (x)dx 

/j 6 S 2 (X)dX 

;; 7 s,(x)dx 

.4 2 

; 3 6 s 3 (x)dx 

/• 7 S 3 (X)dX 
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A’ 

A 23 


A' = 
33 


/^ 7 s 1 (X)dX 

/• 7 S 2 (X)dX 

/; 7 s,(x)dx 

.4 2 

;^ 7 s 3 (x)dx 

/^ 7 S 3 (X)dX 


and 


S (X) is the sensitivity of the yellow layer 

S 2 (X) is the sensitivity of the magenta layer 

S 3 (X) is the sensitivity of the cyan layer 

-The nine "A primes" (above) are just constants dependent only on the 
film sensitivity. The equations (6) through (8) can be solved for EB, EG and 
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ER in terms of E^, E 2 and E^ by several methods. The solutions to equations 
(6), (7) and (8) for EB, EG and ER in terms of E^, and E^ can be written as 


EB - A U E 1 + A 12 E 2 + A 13 E 3 
EG - A n E 1 + A 22 E 2 + A 23 E 3 
ER = A E + A 32 E 2 + A 33 E 3 


(9) 

( 10 ) 

( 11 ) 


These A's are again only a function of the film sensitivity. The 
numerical values of these A’s are given in Table 2 for two different types 
of film. An example of an equivalent energy spectrum is given in Figure 6 


for a blue object. 




ed energy. 


CT 


ER = . 


FILM TYPE 


Kodak 2443 
Kodak 2448 


A 11 

A 12 

A 13 

A 21 

A 22 

A 23 

A 31 

A 32 

A 33 

1.03 

-.0011 

-.02 

-.16 

1.197 

-.032 

-.18 

-.55 

1.7 

1.05 

-.0256 

-.0108 

-.134 

1.176 

-.0155 

-.0008 

-.184 

1.2 


Table 2 


In the case of color IR film, these equations can be written as: 


EG A 11 E 1 + A 12 E 2 + A 13 E 3 


ER A 2i E i + A 22 E 2 + A 23 E 3 
EI = A 31 E 1 + A 32 E 2 + A 33 E 3 


( 12 ) 

(13) 

(14) 


Thus, when E^, E 2 and E 3 have been found, EB, EG and ER can be calculated 
for color film or EG, ER and EI can be calculated for color IR film. The 


physical significance of these equivalent energies, EB, EG, ER and EI, is 
that they represent numbers proportional to average energies over a limited 
band width very similar to the data acquired by multi-spectral scanners. 

EG and ER .correspond exactly to the ERTS satellite MSS bands 4 and 5. 

EI corresponds to a larger band pass than ERTS MSS band 6. 

The calculation for E^ E 2 and E 3 may be difficult if proper procedures 
are not followed. The critical step in the procedure is the exposure of a 
calibrated step wedge on the film before development. The important thing 
here is to be able to calculate E^, and E 3 for each step of the wedge 
from equation (4). In this case the e(A) is the radiation falling on the 
n th step of the film wedge, the film wedge being the exposed step wedge 
on the film of interest. These values are designated as E? (i = 1,3) 
for the step on the film wedge. 


Film Calibration 

The dyes used in color and color IR film are not stable with time. 

On different rolls of film processing differences will cause different 
amounts of dye to be formed in each of the layers for the same exposures. 





For these reasons the familiar D log E or characteristic curve must be 
used for calibration. 


The two ingredients involved in successful calibration of a multi- 
layered film are (1) the step wedge densities as a function of wavelength 
must be known for each step, and (2) the spectral distribution of the 


light incident on the step wedge must be known. 

n th 

Let D (A) = the density of the n step of the step wedge as a 
s 

I 


function of A 

E Q (A) = the incident energy on the step wedge 

- the exposure value for the i th layer from the n th 
of the step wedge 


step 


^ = /c o (X)10 s W S 1 (X)« (15 ; 

1 /&. (A)dA 

1 

If c (A), D (A) , and S.(A) are known, the exposures due to the light 

Q\/J S ^ 

transmitted through each of the steps of the step wedge for each layer can 
be calculated. 

One of the usual practices in sensitometric analysis is to plot the 


density vs (-) the step wedge densities as a characteristic curve. Let 
us look carefully at the assumption one must make for this to be valid. 

If s 0 (A) is a constant, e 0 and D^(A) do not depend on A, then equation 

(15) reduces to 

= e o 10 (16) 

Taking logs of both sides 

log E 1 } = -D^ + log e 0 " (17) 

If the assumptions of D®(A) and e 0 (A) do not depend on wavelengths are true 

the log of the exposure due to the step wedge can be represented by -D^. 

The term log e Q is only necessary if absolute calibration of the system 
is necessary. For most purposes, a relative characteristic curve is 



sufficient, and only the values of the relative densities of the step wedge 
are necessary. It should be noted that in this approximation, the value 
of log t Q is the same for each layer. 

In many cases the density of the step wedge is a constant over the _ 
relevant wavelengths, but the incident energy on the step wedge has a 
distribution like a black body at 5800°K. This is a reasonable assumption 
since most sensitometers are color corrected to this temperature. In this 


approximation, equation (15) reduces to 

fG 0 Si(X) 


e" = 10~ D s 


•>l 5 (e hc/XKT-1) 
/S ± (X) dX 


dX 


(18) 


.n 


E? = 10"°s G. 
1 x 


K s 111 <5 14 (19) 

where G. represents the three different values of ratio of the integrals 

in equation (18). The log of the exposure for each layer then reduces to 

log E? = -D n + log G. (20) 

° l s x 


The important thing to note here is that the term log G^ is different for 
each of the layers. If the exposure for each layer is determined for a 
scene response and these exposures are used independently , only the — 
term need be used. But if the suras of these exposures are used, to 
calculate an equivalent energy spectrum for example, or ratios of the 
exposures are needed, the log G_^ term must be included. 

In the cases that neither d"(X) or e 0 (X) are constant over the relevant 
wavelengths, equation (15) can be used throughout the film analysxs. 


Measurements can be made on each step of the film wedge and dye contents 


constants can be calculated. Three characteristic curves can then be 
plotted; K? T s are plotted against the log of the E^’s as calculated in 
equation (15) . Three different curves (one for each layer) will have 
been plotted. From these curves, for any selected point on the film, values 



of E^, E 2 and can be found for calculated values of K^, K 2 and K^. EB, 
EG and ER or EG, ER and El can be calculated from E^, and for color 
or color IR film, respectively from equations (9) to (14) « 

A necessary ingredient for the calculation of the dye constants K^, _ 
K 2 and K 3 is the knowledge of the spectral densities of the dyes in the 
film. If the shapes of these curves are known and do not change with time, 
the procedures described so far are all that is necessary to calibrate 

t 

multi-layered films. If the height of these curves change with time, the 
step wedge calibration will adequately deal with this problem. But if 
the relative shape of this curve changes with time, or if the spectral 
densities are not known, additional calibrations must be made. 

The easiest way to calibrate the spectral density curves for a film 
is to expose three different step wedges. The best way to expose these 
step wedges is through narrow band interference filters. For color film 
the filters would be centered -at 400 nm, 550 run, and 650 nm.; for color IR 
film the filters would be centered at 550 nm, 650 nm, and 750 nm. 

In the case of color film, the film wedge exposed at 650 nm will 
give the spectral density of the cyan layer. The film wedge exposed 
at 550 nm will contain information on the spectral density of both the 
magenta and cyan layers. If the cyan spectral density has been determined 
then the magenta spectral density can be determined from combination 
density. The film wedge exposed at 400 nm will give the spectral density 
of the yellow dye. A similar argument can be made for color IR film. 

If this technique is employed, the exposures used to plot the 
characteristic curves must be calculated by equation (15) • For color 
film this reduces to: 

E n /e o (400)10 _D s C400) F 1 S 1 (400)dA (2i; 


/Sj^JdA 



( 22 ) 


,,0 _ e (400)l(f D s (400) S, (400) 
/ S 1 (X)dX 

“ e o (550)10’’ D s (550) S 2 (550) 
/S 2 (A)dX 

“ e o (650)10 " D s (650) S 3 (650) 
/S 3 (A)dA 


(23) 


(24) 


The characteristic curve for is determined from the film wedge exposed 
at 400 nm; the curve for K ? is determined from wedge exposed at 550 nm; 
and the curve for is determined from the wedge exposed at 650 nm. The 
curve for K 3 could be determined from the wedge exposed at 550 nm also. 

In that case, would be: 

F n _ e o (550)KfV 550) S 3 (550) (25 

3 /S 3 (X)dX 

An alternate method of arriving at the same result is to expose one 
step wedge with a sensitometer and also one step of the wedge through the 
three different narrow band filters. Using the step exposed through the 
filters the spectral densities of the dyes can be determined. The charac- 
teristic curves are then plotted from the complete step wedge exposed by 
the sensitometer. 


Summary 

We have presented a method of quantitatively analyzing multi-layered 
films. This type of analysis uses the physical characteristics of the film 
to relate the film’s spectral response to the spectral distribution of the 
incoming radiation. A method has been presented to determine three dye 
content constants at each point on the film that can be easily related 
to the incoming radiation. The equivalent energies are simply a linear 
transformation of these dye content constants which are more closely 
related physically to the incoming radiation. These equivalent energies 



can be thought of as the output from a three channel multi-spectral 
scanner. In this regard this method is well suited for comparison to 
ERTS satellite data. 


REFERENCES 

1. Kodak Data for Aerial Photography, M-29 , Eastman Kodak Company, 

1971, 62p . 

i 

2. Ibid . 

3. Ibid . 

4. Manual of Color Aerial Photography , John T. Smith, Jr., (Ed.), 

American Society of Photogrammetry , 1968, 308p, 

5. Brauer, Nohel, Schneider, Linear Mathematics , W. A. Benjamin, Inc. 


New York, 1970. 



